Dynamic effects on the loss of control in template-directed nucleation 
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Full nucleation control for deposited functional molecules on pre-patterned surfaces is of major 
technological relevance. To understand the nucleation behavior we combine the numerical solution 
for the evolution of the adatom concentration with standard nucleation theory. From the qualitative 
change in nucleation behavior upon variation of the pattern spacing and coverage we show why 
the quality of nucleation control can vary significantly in different parameter regimes. In some 
limits analytical expressions can be formulated for the nucleation control. Our analysis provides a 
theoretical explanation for previous experimental observations [Wang et al, PRL 98, 225504 (2007)]. 
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The controlled fabrication of nano-structures is a tech- 
nologically highly demanding and interesting topic. In 
this context techniques based on lithography and etch- 
ing have become very sophisticated over the past decades. 
However, with the advent of functional organic molecules 
in micro-electronics, one faces the problem that this tech- 
nique is often not applicable to this class of substances. 
At this point the concept of template-directed growth be- 
comes interesting: in this technique the areas of desired 
accretion are marked in some fashion to favor the adsorp- 
tion of the functional substance, which is subsequently 
deposited. There exists a number of methods to generate 
these preferred sites; see, e.g., Refs. P, 0, H, 0, H, H, H Hi • 

From a technical point of view one is interested in ex- 
clusive accretion of the deposited substance at the pre- 
defined sites and no additional nucleation. The prede- 
fined sites display a regular quadratic arrangement (size 
of unit-cell: p 2 ; see inset of Fig. [TJ. To provide a quanti- 
tative means of evaluation, the nucleation control can be 
expressed by 
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with R representing the number of additional nucleated 
clusters per unit-cell. Additionally, to be able to relate 
theory and experiment the length scales are renormalizcd 
by 




to p := — 



P 
A ' 



(2) 



where (N/A) n represents the overall island density of an 
unpatterned substrate. For p 3> 1 the influence of the 
predefined sites on the nucleation processes can be ne- 
glected so that the nucleation properties just follow the 
standard nucleation theory; see, e.g.,0, E3, EH, E2, Ell- 
in contrast, for dense patterns with p <C 1 all nucleation 
occurs at the predefined sites (xn ~ 1) 0, E3|- 



The transition regime p « 1 is of particular interest. 
First, from an experimental/technological view one needs 
to understand how far the range of complete nucleation 
control extends to large pattern sizes. Second, from a 
theoretical perspective new phenomena come into play as 
compared to standard nucleation theory because of the 
interplay of the nucleation tendency and the drainage of 
the adatoms by the predefined sinks. The complexity of 
this regime is also highlighted in Fig. [D where the results 
of previous computer simulations [a 1 141 ] and experiments 
are displayed [8( . Whereas the simulations display nucle- 
ation control in the regime p > 1 the experiments show a 
loss of nucleation control at pattern spacings significantly 
below A (i.e. p < 1). What is the interpretation of these 
differences? 
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FIG. 1: Comparison of experimental (x) and simulated (O) 
nucleation control in for different pattern spacings p from 
Ref. [U]. While the simulated curve shows a retention of nu- 
cleation control for p > 1, the experimental curve exhibits an 
early loss at p significantly smaller than unity. In the inset 
the parameters of the model are defined. 

In this paper we present the theoretical basis of nu- 
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cleation control and in particular reveal the influence of 
the experimental deposition time scale. For a clear iden- 
tification of the physical mechanisms we employ a mini- 
mum model which considers particle deposition at flux F 
and subsequent surface diffusion with diffusion constant 
D. Adatoms are either included by the predefined sites, 
acting as perfect sinks, or, if the distance between the 
sites p is too large or F/D is too small, start to form 
additional nuclei (i?) aside of these locations. Since we 
are interested in the key mechanisms, we neglect possible 
desorption processes. 

To evaluate the temporal evolution of the nucleation 
behavior on a pre-patterned substrate a simulation tech- 
nique similar to the level-set methods of Refs. fljl liL 17 . 
flil | was applied. This algorithm is far more efficient as 
compared to the previous atomistic Monte Carlo simula- 
tions [1, • It is based on the numerical solution of the 
adatom diffusion equation 



<9c(x) 
dt 



L>V 2 c(x) + F 



by discretization in time and space. The generation of 
new nuclei follows the concepts of standard nucleation 
theory 0, El E2- Thc local 

rate and hence also thc 



chance of nucleation Sjj on the grid point j at time t 
with the next time step At is given by 

Sj,t = fc(c 3 - t /c ) r+1 AxA 2 /Ai, 



(4) 



where i* + 1 represents the number of adatoms in the 
first stable two-dimensional nucleus and k and cq are pa- 
rameters which characterize the nucleation probability. 
The parameters [l9| were adjusted such that we obtained 
a typical distance of about 50 grid points between nu- 
clei in the steady state regime of the unpatterned case. 
For the definition of a length scale in our simulations we 
have introduced the lengthscale u, where u/2 describes 
the distance of two adjacent grid points. Eq. (J3]) is sup- 
plemented by the boundary condition c s i„k = at the 
positions of nuclei. 

Commonly (see e.g. [Ill [lil, [l?]] ) , a new nucleus is gen- 
erated when Q(t), defined via Q = J2j t &j>t excee ds an 
integer value. However, the deterministic nature of the 
first nucleation event introduces a bias during a critical 
phase of this simulation [2(J. Thus, we apply the fol- 
lowing Monte Carlo like algorithm: at every time step, 
St = J2i $i,t is spatially integrated and in combination 
with sufficiently small At provides St At -C 1. This 
chance is then compared to a random number < r < 1 
resulting in nucleation for r < St At. In a successful event 
the new nucleus is placed stochastically according to the 
individual Si t f 

A simulation run consists of one quadratic unit-cell of 
the surface pattern with length p as shown in the inset 
of Fig. [U using periodic boundary conditions. We have 
checked that the xjy vs. p-plot does not display any finite- 
size effects. The finite-size effects, reported in literature 



2l| . occur in a regime where in the present case ijy = 1. 
Each data point results from averages over 200 indepen- 
dent simulations. All systems are subject to the same 
flux F and a value of i* 



4 [14|. Similar to Refs. [22j,|2 



we represent the predefined sites as a single grid point 
on our discretized simulation lattice, acting as a perfect 
sink. The same holds for newly generated nuclei. In this 
way a predefined or generated nucleus possesses an effec- 
tive area u 2 /4. As we are not interested in the temporal 
evolution of the individual nuclei, but in the averaged 
properties of the surface, the sinks remain to be confined 
to their initial position and do not expand any further. 
Experimentally, the area Q occupied by an adatom is, 
of course, much smaller than that of the pattern. Con- 
sequently, the ratio fi/u 2 is orders of magnitude smaller 
than unity (0(1O -4 )) and the resulting values for the cov- 
erage 0, which are given below in units of fi/u 2 , hence 
correspond to sub-monolayer growth. 

For good nucleation control the nucleus density of the 
(3) pre-patterned surface (N/A) pt has to be lower than that 
of the unpatterned counter part (N/A) n . In what follows 
we compare the temporal evolution of these densities for 
several patterned systems of different p complemented 
by an unpatterned reference system; see Fig. [21 The un- 
patterned system (solid line) possesses the characteris- 
tic shape expected from standard nucleation theory [Io| : 
The beginning shows a transient time regime of rapid 
nucleation, dominated by local supersaturation. How- 
ever, as the surface becomes saturated with nuclei, the 
overall concentration decreases and the system enters the 
regime of steady state nucleation. Mechanistically, the 
comparatively slow nucleation rate of this regime can be 
attributed to statistical fluctuations in the adatom con- 
centration. 
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FIG. 2: Evolution of the nucleus densities N/A with growing 
time. Shown are surfaces of different pre-patterned densities 
(with p/u indicated by the numbers), as well as the unpat- 
terned reference system (solid line). 



With respect to the patterned surfaces Fig. fallows the 
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following observations: (a) For comparably large p (see, 
e.g, p/u = 90) the behavior of the patterned systems re- 
sembles that of the unpatterned counter part, but as by 
construction there already exists one dot at the beginning 
of the simulation, (b) the systems start off at a density of 
(N/A) pt {t = 0) > (N/A) u (t = 0) = 0. As growth com- 
mences, however, the permanent adatom depletion by 
the pattern leads to an effectively smaller nuclcation rate 
such that (c) the patterned surfaces approach the unpat- 
terned long time behavior from below, (d) Furthermore, 
as p becomes smaller (p/u = 50, 40) the s-shaped part of 
the curve is stretched in time corresponding to a slower 
transition through the rapid nucleation regime. This ef- 
fect increases to the point that (e) at sufficiently low 
pattern spacing, e.g. p/u = 32.5, the s-shape becomes 
unrecognizable and (N/A) pt remains significantly below 
(N/A) u within the considered time scale, (f) Finally, p 
can be decreased to the point where (N/A) pt (t = 0) be- 
comes so large that the unpatterned long time behavior 
is approached from above. 

Returning to the question of nucleation control, Eq. ([I]) 
can, with the help of Eq. be rewritten as 




(5) 



p 2 and show re- 



A surface can hence only lie above xn 
tained nuclcation control if (N/A) u /(N/A) pt > E Specif- 
ically this means that because of (b) the nucleation con- 
trol at the beginning of the experiment is always very 
low, but due to (c) and (d) this changes with increasing 
time for all patterns up to (f). 

The corresponding nucleation control can be found 
in Fig. [3] at low coverage (9 = 60f2/u 2 ) the unpat- 
terned surface does not yet show significant nucleation 
and hence the vast majority of patterned systems exhibits 
a higher nucleus density, which in turn results in an early 
loss of nucleation control, indicated by a curve in the 
lower left segment of the figure. As time increases, how- 
ever, and the regime of rapid nucleation begins around 
9 = 70£l/u 2 , the unpatterned nucleation density com- 
mences to exceed its patterned counterparts leading to 
an increasing number of systems under nucleation control 
and hence to a crossing of the in = p -2 diagonal closer to 
p = 1. As the unpatterned surface reaches the beginning 
of the steady state nucleation regime, nucleation control 
becomes most pronounced which can be seen from the 
curve corresponding to 9 = 100il/u 2 . Past this point 
the slow convergence of all patterned nucleus densities to 
that of the unpatterned surface results again in a decrease 
of nucleation control. All this effects together mean that 
at fixed p one has a non- monotonous dependence of xjv 
on O. Note that upon changing the coverage the onset 
of nucleation control can vary as much as a factor of two 
with respect to the critical pattern density (relative to 
the respective value of A). 



.f'N 



0.1 - 



□ 


□ 60 





O 100 


A 


A 350 



•A 



''A 



p 



FIG. 3: Nucleation control in with increasing coverage 
9/(fi/u 2 ) = 60(D), 100 (O) and 350(A). At relatively low 
coverage the patterns display an early loss of nucleation con- 
trol (in)- However, as the unpatterned reference system 
passes into the steady state nucleation regime, nucleation con- 
trol becomes most pronounced. With further increasing time 
the curves converge to the dotted limit. The low-coverage 
limit :rjv = l/(l+p 2 )is indicated by the broken line. 



Two limits can be treated more quantitatively. First, 
in the limit of low coverage and large p (Op 2 = const) 
the growth of additional nuclei is not influenced by the 
gold dot. Thus one has (N/A) pt = (N/A) u + p~ 2 , which 
yields xn = 1/(1+ p 2 )- This limit is included in Fig. [3] 

Second, further analysis is possible when the first nu- 
cleus is formed from a stationary concentration profile. 
As shown in Ref. 14| for fixed parameters one can find 
a critical time scale t c such that for t > t c the station- 
ary concentration profile is reached. For this condition 
to be applicable it is important that nucleation can be 
neglected for times smaller than t or, equivalently, R < 1 
at time t. The time scale t c scales like p 2 [3]. In this 
scenario one can derive the relation 



Roc&p 



(6) 



where / can be calculated from knowledge of the sta- 
tionary concentration field, which was done numerically 
in Ref. Q. 

In Figured] we have replotted Xn{p) as R(p). In this 
representation one can clearly see that Eq[5] is recovered 
for small values of p. This is to be expected because 
for a given time scale of the experiment the stationarity 
condition is always fulfilled in the limit of small p and 
correspondingly small t c . Furthermore the linearity with 
respect to 9 is reflected in the linear increase of (N/A) for 
small p in Fig. [2] Interestingly, the applicability range of 
this power law strongly depends on the coverage: while 
at 9 = 60£!/u 2 the relation only holds for values less 
than R w 0.05, for larger coverage the power law nearly 
holds up to the theoretical limit of R = 1. This can 
be easily understood from analyzing R at the time scale 
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t c , introduced above. Joining together Eq. ([6]) and the 
relation t c oc p 2 one obtains R{t c ) oc t^J 2 . Indeed, one 
observes a dramatic increase of the range of applicability 
of Eq. © upon increasing the coverage. Of course, at 
the latest the power law scaling has to break down when 
R approaches 1. 
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FIG. 4: Additionally formed nuclei R(p) at increasing cover- 
age e/(fi/u 2 ) = 60 (□), 100 (O) and 350 (a). The theoretical 
expectation for the static case according to Ref. [141 ] (R oc p' 
with / = 13.2) is represented by the dashed lines. The regime 
of retained nucleation control is indicated in grey. 

Based on these considerations we can conclude the fol- 
lowing: the experimentally observed early loss of nucle- 
ation control as depicted in Fig. [T] does not disagree with 
the general understanding of the mechanisms governing 
template-directed growth, and can be attributed to dy- 
namic effects. In particular, this effect is caused by a dif- 
ferent temporal evolution of the nucleus density on the 
patterned and unpatterned substrates: unpatterned sub- 
strates or patterns of very large p show a distinct separa- 
tion between the transient and the steady state regime, 
which mechanistically corresponds to nucleation by su- 
persaturation or by statistical aggregation, respectively. 
With decreasing pattern spacing, however, the transient 
regime becomes dilated until it is phenomenologically as 
well as mechanistically indiscernible from steady state 
nucleation. In the commonly used representation, which 
is normalized by an unpatterned surface, this leads to the 
effect that initially the majority of patterned systems ex- 
hibits an early loss of nucleation control. This changes, 
however, as the unpatterned surface approaches the end 
of the transient regime which leads to maximal nucleation 
control close to the beginning of steady state nucleation. 
Thus, a consistent picture of the influence of the distance 
of the patterns as well as the experimental time scale on 
nucleation control can be formulated. This is of partic- 
ular relevance for the understanding of the formation of 
nanostructure on surfaces. Evidently, also more complex 
situations such as non-circular patterns can be analyzed 



within the present model approach. 
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